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ABSTRACT 

An upwind 3-D finite volume Navier-Stokes code is modified to facilitate modeling of 
complex geometries and flow fields presented by proposed National Aero-Space Plane concepts. 
Code enhancements include an equilibrium air model, a generalized equilibrium gas model, and 
several schemes to simplify treatment of complex geometric configurations. The code is also 
restructured for inclusion of an arbitrary number of independent and dependent variables. This 
latter capability is intended for eventual use to incorporate nonequilibrium/chemistry gas models, 
more sophisticated turbulence and transition models, or other physical phenomena which will 
require inclusion of additional variables and/or governing equations. Comparisons of computed 
results with experimental data and with results obtained using other methods are presented for code 
validation purposes. Good correlation is obtained for all of the test cases considered, indicating the 
success of the current effort. This work was conducted at the NASA Langley Research Center, 
during participation in the NASA/Industry Fellowship Program for the National Aero-Space Plane. 


*This work was conducted at the NASA Langley Research Center, during participation in the 
NASA/Industry Fellowship Program for the National Aero-Space Plane 
t Engineering Specialist 
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INTRODUCTION 


The National Aero-Space Plane (NASP) program has highlighted the need for development 
of advanced computational fluid dynamics methodology. The success of the program, unlike that 
for any previous aircraft, depends upon the availability of the state-of-the-art in flow simulation 
and prediction. Advances in flow discretization techniques, solution algorithms, equilibrium and 
nonequilibrium/chemistry gas models, and turbulence and transition models must be incorporated 
into methodology capable of treating the complex geometries and flow fields presented by 
proposed NASP concepts. 

The NASA/Industry Fellowship Program provided this author with an opportunity to assist 
in the development of one such method. The basic CFL3D code, an advanced thin-layer Navier- 
Stokes flow solver which is relatively easy to use and which features the flexibility required to treat 
complex flows, was modified during this effort to incorporate equilibrium air and generalized 
equilibrium gas models, and to further enhance its geometric modeling capabilities. At the same 
time, the code was restructured to facilitate future computations incorporating an arbitrary number 
of independent and dependent variables. This latter capability is intended for eventual use to 
incorporate nonequilibrium/chemistry gas models, more sophisticated turbulence and transition 
models, or other physical phenomena which will require inclusion of additional variables and/or 
governing equations. 


NOMENCLATURE 


a 

speed of sound 

e 

internal energy per unit mass 

L 

reference length 

M 

Mach number 

P 

pressure 

Pr 

Prandtl number 

Re 

Reynolds number 

T 

temperature 

u,v,w 

Cartesian velocity components 

x,y,z 

Cartesian spatial coordinates 
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y + wall unit. Ay ■ (p w Ohv I Pw) ,/2 

y specific heat ratio 

y "equivalent" specific heat ratio, 1 + pi pe 

r "equivalent" specific heat ratio, a 2 pi p 

p mass per unit volume (density) 

k thermal conductivity 

p viscosity 

to vorticity 

Subscripts 


oo freestream 

w wall 


DESCRIPTION OF METHOD 

The computer program to be described is derived from the April 1988 release of CFL3D 
(Version 1.0), a method which is well documented in the open literature [1,2,3, 4]. A brief outline 
of CFL3D methodology is given below, followed by a discussion of enhancements and features 
incorporated in the present code. 

Overview of Basic CFL3D Methodology 

The governing flow equations are the three-dimensional, time-dependent, conservation law 
form of the compressible Euler or thin-layer Reynolds-averaged Navier-Stokes equations, 
expressed in generalized coordinates. An upwind-biased approach with up to third order accuracy 
is used to evaluate the inviscid fluxes at the cell interfaces, as described below. A spatially-split, 
three-factor approximate factorization algorithm and Euler implicit time integration/linearization is 
used to advance the solution (cell-averaged flow properties) in time [5]. 

Inviscid flux interface values are obtained using a MUSCL interpolation scheme [6], 
coupled to either the flux difference splitting (FDS) scheme of Roe [7,8] or the flux vector splitting 
(FVS) scheme of Van Leer [5,9]. Flux splittings are based on a one-dimensional Riemann 
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problem, and are subsequently modified to treat multi-dimensional flows. Overall, these 
approaches provide an upwind-biasing in the flux interface evaluation. They also introduce an 
amount of dissipation which is consistent with the discretization of the governing flow equations, 
and which is required to stabilize the solution procedure. The so-called smooth and min-mod flux 
gradient limiters are optionally employed, to minimize the adverse effects of large flow gradients 
and discontinuities (such as shock waves). 

At each time step, the FDS/FVS approaches lead to a series of 5-by-5 block tridiagonal 
matrix inversions, for each of the spatial directions. Additional approximations may also be made 
in the FDS scheme so as to diagonalize the solution matrices [3]. This leads to a series of scalar 
tridiagonal matrix inversions, and an attendant reduction in execution time. 

Viscous and heat flux interface values are obtained using central finite-difference formulae. 
The laminar thin-layer Navier-Stokes terms may be included in all three directions. A Baldwin- 
Lomax algebraic turbulence model [10] is also employed. The effects of turbulence may be 
included in one direction, or in two directions via a distance-weighted two-wall comer model for 
the turbulent eddy viscosity. 

A zonal grid structure facilitates modeling of complex geometries and/or flow fields. 
Explicit treatment of grid boundaries further simplifies this task, since boundary condition subrou- 
tines are easily modified for specific or unusual cases. Provisions are included for treatment of 
blocked grids, longitudinally-patched grids [1 1], and dynamic moving grids. A variety of cell-cen- 
ter or cell-interface type boundary conditions may also be specified at grid boundaries ... 
freestream flow, extrapolation from the interior (supersonic outflow), subsonic characteristic in- 
flow/outflow (based on one-dimensional Riemann invariants [12]), inviscid wall (flow tangency), 
viscous wall (adiabatic or fixed wall temperature), and an assortment of 
symmetry/periodicity/singular-axis/wake-continuation type boundary conditions. 

Several schemes are available to reduce overall execution time, particularly for computing 
steady flows. Local-time-stepping and multigrid [12,13,14] techniques accelerate code 
convergence. Mesh sequencing is a technique whereby solutions obtained on coarser grids are 
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used to initialize flow field data on successively finer grids, until finally a solution is obtained on 
the desired input grid. Mesh embedding is a technique whereby enhanced solution accuracy is 
obtained by locating even finer grids in particular regions of interest. Both mesh sequencing and 
mesh embedding reduce the computational effort expended to achieve a given level of solution 
accuracy, and their use is facilitated by automated grid generation and flow field interpolation 
routines. 

In February 1989, while the present code was still under development, an updated version 
of CFL3D (Version 1.1) became available. The enhanced capabilities of the updated version of 
CFL3D were subsequently incorporated in the present code, including an improved treatment for 
longitudinally-patched grids [15], more generalized boundary conditions, and mesh sequencing for 
two-dimensional flows. The only CFL3D enhancement not found in the present code is an 
alternate two-factor approximate factorization algorithm [12], 

Equilibrium Air and Generalized Equilibrium Gas Models 

The methodology described above assumes a perfect gas model for the thermodynamic and 
transport properties of the fluid. Versions 1.0 and 1.1 of CFL3D further assume air to be the 
working fluid. These restrictions do not apply to the present code, which features more general 
equilibrium gas capabilities. 

The flux-splitting schemes of Roe and Van Leer are extended in the present code to treat 
real gases, using techniques developed by Grossman and Walters [16], The perfect gas 
relationships are replaced by equilibrium gas relationships, usually in the form of curve fits. The 
specific heat ratio y employed in the flux-splitting schemes is then replaced by the "equivalent" 
values y= 1 + p / p e and r = a 2 p/ p. 

Two equilibrium gas models for thermodynamic properties (p, p , c, a and T) are 
incorporated into the present code. The first, due to Srinivasan and Tannehill [17], consists of 
curve fits for equilibrium air, and executes in scalar mode. The second, due to Liu and Vinokur 
[18], is a generalized equilibrium gas model, uses bicubic spline interpolation (based on an 


5 


auxiliary interpolation coefficient data file), and executes in vector mode. An interpolation 
coefficient file for equilibrium air obtained from Liu was augmented by this author in order to use 
the approach at lower temperatures normally considered to be in the perfect gas regime. 

The equilibrium gas model for transport properties (p, k, and Pr) is the equilibrium air 
curve fits due to Srinivasan and Tannehill [19]. Versions which execute in vector mode were 
developed by this author, after discovering that more execution time was used for computing 
transport properties than for computing thermodynamic properties. 

The Liu and Vinokur thermodynamic property model is not restricted to equilibrium air, 
since auxiliary interpolation coefficient data files for other equilibrium gases could be constructed. 
A similar approach for the transport properties is hopefully underdevelopment, and when available 
can be incorporated into the present code as well. 

Relative to perfect gas computations, the original (scalar) Srinivasan and Tannehill 
thermodynamic and transport property gas models result in roughly a 125% increase in execution 
time. Using the vectorized Srinivasan and Tannehill transport property model results in only about 
a 50% increase in execution time, while using the vectorized Srinivasan and Tannehill transport 
property model and the vectorized Liu and Vinokur thermodynamic property model results in only 
about a 20% increase in execution time. Of course, these numbers are approximate, and reflect 
average values obtained for a variety of test cases. 

First Steps towards a More Generalized Flow Solver 

In conjunction with the equilibrium gas flux-splitting capability, the present code was 

enhanced so as to permit an arbitrary number of independent and dependent variables to be stored 
in the ^-vector. 

The basic CFL3D code stores only the five independent variables p, u, v, w, and p in the 
^-vector. In the present code, an arbitrary number of independent variables (i.e., the number of 
conserved variables or governing equations), Iqcv, and an arbitrary number of dependent variables 
(e.g., yand T), Iqdv, may be stored in the ^-vector. The values of Iqcv and Iqdv need be set only 
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once, in a parameter statement in the main program. These values, and the total number of 
variables in the q- vector, Iqt = Iqcv + Iqdv, are subsequently passed to the required subroutines as 
arguments and/or through common blocks, for appropriate dimensioning of arrays and indexing of 
do loops. 

This coding structure is a first step towards a more generalized flow solver which might 
incorporate nonequilibrium/chemistry gas effects, more sophisticated turbulence and transition 
modeling, or other physical phenomena which will require inclusion of additional variables and/or 
governing equations. Additional work will be required before this goal is achieved. For example, 
matrix inversion logic is currently fixed to treat 5 governing equations, and increasing the value of 
Iqt may result in overlap or overflow of flux routine scratch arrays. Nevertheless, the majority of 
the present code should not require further modification in order to incorporate more generalized 
flow models. 


Other Code Enhancements 

Two important features of the present code were developed to enhance user friendliness. 
First, the path and name of all auxiliary data files (currently as many as 11) are specified via the 
standard unit 5 input data, rather than in the FORTRAN coding itself, to avoid code modification 
and recompilation. Second, user specified scale factors for length, mass, and temperature permit 
the use of arbitrary dimensions (e.g., metric or English) for the input and output data. 

Other features of the present code offer enhanced capabilities. Most significant of these is a 
very generalized grid blocking boundary condition capability (developed by George Switzer, 
Analytical Services and Materials, Inc.). Also noteworthy is a "jagged" boundary condition 
algorithm (developed by Mark Eppard, Analytical Services and Materials, Inc.) which permits 
treatment of surface edges that are skewed with respect to, or cut across, grid lines. A new flux 
interface averaging procedure, developed at NASA Langley, may enhance convergence for cold 
wall cases. Since the jagged boundary condition and flux interface averaging capabilities are not 
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yet fully generalized, they have been commented out in the FORTRAN coding (lines start with the 
characters "ebsr"), and should be activated only by knowledgeable users. 

RESULTS 

Results computed for several test cases are presented in order to evaluate the present code's 
capabilities. For each test case, calculations are compared to results obtained using other methods, 
or to experimental data, and previous comparisons by other investigators are cited. 

Computations were obtained using the perfect gas model, the Srinivasan and Tannehill 
equilibrium air model, and the Liu and Vinokur generalized equilibrium gas model with the 
augmented auxiliary interpolation coefficient data file for equilibrium air. Since the two 
equilibrium air models gave essentially identical results for all of the test cases, only those obtained 
with the Liu and Vinokur model are presented herein. 

Unless otherwise noted, all results were computed using FDS, third order upwind-biased 
spatial accuracy, min-mod flux limiter, and the 5-by-5 block tridiagonal matrix inversion 
algorithm. The majority of the computations were made for laminar flow, and included thin-layer 
terms in the ^-direction (normal to the body surface) only. A fixed wall temperature was specified 
for use in all viscous wall boundary conditions. Local time stepping was used to accelerate 
convergence to steady state. 

Supersonic Laminar Flat Plate Boundary Layer 

The first test case consists of supersonic laminar flow over a flat plate (this is also one of 
the test cases studied in [20]). The flow conditions are: 

Moo = 2.0 

Reoo/L = 1.65-10 6 /m 

Too = 221.6 °K 

T w = 221.6 °K 
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A grid consisting of 51 grid points in the streamwise direction and 100 grid points normal to the 
surface was employed. Average grid spacing normal to the surface was 0.43- 10" 4 m, producing 
an average y + of 1.33. The residual was reduced approximately 4.5 orders of magnitude over 
4000 time steps. NASA Cray-YMP (Reynolds) execution times required for the perfect gas, 
Srinivasan and Tannehill, and Liu and Vinokur gas models were 4.5- 10~ 5 , 7.2- 10~ 5 , and 
5.4- 10 -5 cpu-seconds per mesh-cell-point per time-step-iteration, respectively. 

Computed supersonic laminar flat plate boundary layer results are compared to predictions 
made using a conventional boundary layer calculation [21] (boundary layer calculations supplied 
by Douglas Dilley, Analytical Services and Materials, Inc.). Velocity and temperature profiles at 
an axial location x = \m are presented in Fig. 1. Axial distributions of heat transfer and skin 
friction are presented in Fig. 2. All of the present results show excellent correlation with the 
boundary layer predictions. As expected, equilibrium gas effects are not significant for this 
relatively low temperature flow. 

Hypersonic Laminar Flat Plate Boundary Layer 

The second test case consists of hypersonic laminar flow over a flat plate (this is also one 
of the test cases studied in [22]). The flow conditions are: 

A/co = 20.0 

ReootL = 2.010 5 / m 

Too = 100.0 °K 

T w = 1000.0 °K 

A grid consisting of 64 grid points in the streamwise direction and 64 grid points normal to the 
surface was employed. Average grid spacing normal to the surface was 0.1- 10 -3 m, producing an 
average y + of 1 .08. The residual was reduced approximately 5 orders of magnitude over 4500 
time steps. Reynolds execution times required for the three gas models were 4.4- 10~ 5 , 8.810 -5 , 
and 5.5 10 -5 cpu-seconds per mesh-cell-point per time-step-iteration, respectively. 
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Computed hypersonic laminar flat plate boundary layer results are compared to predictions 
made using CFL3DE, an extension of the CFL3D method by other investigators [23] which also 
incorporates equilibrium air effects (CFL3DE calculations supplied by Douglas Dilley, Analytical 
Services and Materials, Inc.). Velocity and temperature profiles at an axial location x = Im are 
presented in Fig. 3. Axial distributions of heat transfer, skin friction, and pressure are presented in 
Fig. 4. The present results show excellent correlation with the CFL3DE calculations. Equilibrium 
gas effects are significant, particularly for the temperature profile predictions. 

High Speed Inlet 

The third test case is the high speed flow through an inlet (this is also one of the test cases 
studied in [16]). The flow conditions are: 

Moo = 5.0 

Reoo/L = 4.9410 6 / m 

Too = 3573.0 °K 

The inlet features a 10° compression, followed downstream by a 10° expansion. Inviscid 
computations were obtained, to permit comparison with the exact perfect gas and equilibrium air 
solutions. A grid consisting of 201 grid points in the streamwise direction and 51 grid points 
normal to the surface was employed. The residual was reduced approximately 3 orders of 
magnitude over 3000 times steps. NASA Cray-2 (Navier) execution times required for the three 
gas models were 1.310 -4 , 1.610 -4 , and 1.3- 10 -4 cpu-seconds per mesh-cell-point per time- 
step-iteration, respectively. 

Computed high speed inlet results are compared to the exact inviscid solutions. Inlet-wall 
density, pressure, and temperature distributions appear in Fig. 5. The agreement is good, except 
for the temperature level aft of the expansion, which is overpredicted. The same effect is seen in 
[16]. No attempt was made to try to eliminate the post-shock oscillation evident in the present 
predictions, which nonetheless indicate the proper perfect gas/equilibrium air trends. 
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Calculations were also made using FVS. The residual was reduced 3.5 orders of 
magnitude over 3000 time steps. Navier execution times required for the three gas models were 
8 .MO- 5 , 1.2- 10~ 4 , and 8.3- 10~ 5 cpu-seconds per mesh-cell-point per time-step-iteration, 

respectively. The results are shown in Fig. 6, and are similar to those obtained using FDS. 

Bent Nose Biconic 

The fourth test case is high speed laminar flow past a bent nose biconic (one of the test 

cases studied in [24]). The flow conditions are: 

Moo = 9.86 

ReoJL = 1.842- 10 6 /m 

Too = 49.75 °K 

T w = 300.0 °K 

As shown schematically in Fig. 7, a total of 85 grid points in the streamwise direction, 45 grid 
points normal to the surface, and 23 grid points circumferentially was used to model one-half of 
the configuration, with symmetry imposed across the x-z plane. Average grid spacing normal to 
the surface was 0.5- 10~ 5 m, producing an average y + of 0.23. To avoid difficulties sometimes 
encountered using FDS to compute blunt nose flow fields, FVS was employed. The mesh 
sequencing capability was also used, to minimize overall execution time. The residual was reduced 
approximately 4.5 orders of magnitude over 4300 time steps. NASA Cray-2 (Voyager) execution 
times required for the three gas models were 7.0T0 - ^, 9.9- 10~^, and 8.3-10 ^ cpu-seconds per 
mesh-cell-point per time-step-iteration, respectively. 

Computed bent nose biconic surface heat transfer rates are compared to experimental data 
[25] in Fig. 8. The present results show good correlation with the data. Equilibrium gas effects 
are less significant than expected for this high speed flow. 


11 



Flared Cone (Laminar) 


The fifth test case is that of high speed laminar flow past a flared cone (one of the test cases 
studied in [24]). The flow conditions are: 

Moo = 16.93 

Reoo/L = 1.976 10 5 /ft 

Too = 83.73 °R 

T w = 530.0 °R 

As shown schematically in Fig. 9, a total of 97 grid points in the streamwise direction, 45 grid 
points normal to the surface, and 19 grid points circumferentially was used to model one-half of 
the configuration, with symmetry imposed across the x-z plane. Average grid spacing normal to 
the surface was 0.241(H/r, producing an average y + of 0.09. Employing the mesh sequencing 
capability, the residual was reduced approximately 3.5 orders of magnitude over 3100 time steps. 
Navier execution times required for the three gas models were 1.2-10 -4 , 1.610- 4 , and 1.5- 10~ 4 
cpu-seconds per mesh-cell-point per time-step-iteration, respectively. 

Computed flared cone surface heat transfer, skin friction, and pressure distributions for 
laminar flow are compared to experimental data [26] in Fig. 10. The present results show good 
correlation with the data. Equilibrium gas effects are less significant than expected for this high 
speed flow. 


Flared Cone (Turbulent) 

The sixth test case considered is high speed turbulent flow past a flared cone. The flow 
conditions are: 

Moo = 7.85 

Reoo/L = 4.697- 10 6 /ft 

Too = 130.2 °R 

T w = 530.0 °R 

The grid, shown schematically in Fig. 11, is similar to that for the laminar case. Average grid 
spacing normal to the surface was 0.83- 10~ 5 /t, producing an average y + of 0.51. The turbulence 
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model was employed in the ^-direction (normal to the body surface) only. As for the laminar case, 
mesh sequencing was employed, and the residual was reduced approximately 3.5 orders of 
magnitude over 3100 time steps. Voyager execution times required for the three gas models were 
6.410 -5 , 8.710 -5 , and 7.210 -5 cpu-seconds per mesh-cell-point per time-step-iteration, 

respectively. 

Computed flared cone surface heat transfer, skin friction, and pressure distributions for 
turbulent flow are compared to experimental data [26] in Fig. 12. The present results show good 
correlation with the data. The figure clearly indicates that appropriate use of the algebraic 
turbulence model can enhance code predictions. 

Laminar Comer Flow 

The seventh test case consists of laminar flow in a comer formed by two intersecting 
wedges (this flow is also studied in [27]). The flow conditions are: 

Moo = 3.0 

Reoo = 2 . 22 - 10 5 

= 105.0 °K 

T w = 294.0 °K 

The flow was computed on a 120 by 120 crossflow plane grid, assuming conical flow in the 
streamwise direction. Average grid spacing normal to the surface was 0.1410 -3 times x, 
producing an average y + of 6.05. Laminar viscous thin layer terms normal to both walls were 
included, in the j- and ^-directions. The residual was reduced approximately 3 orders of 
magnitude over 3600 time steps. Navier execution times required for the three gas models were 
7.810 -5 , 1.210 -4 , and 1.M0 -4 cpu-seconds per mesh-cell-point per time-step-iteration, 

respectively. 

Computed wall pressure distributions for laminar comer flow are compared to experimental 
data [28] in Fig. 13. The present results show good correlation with the data. 
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Turbulent Corner How 


The eighth test case consists of turbulent corner flow (also studied in [27]). The flow 
conditions are: 


Moo 

= 3.0 

Reoo 

= 3.0310 6 

Too 

= 105.0 °K 

T w 

= 294.0 °K 


The flow was again computed on a 120 by 120 crossflow grid, and was assumed to be conical. 
Average grid spacing normal to the surface was 0.1010 -3 times x, producing an average y + of 
5.51. The two-wall comer model was used to simultaneously include turbulence effects normal to 
both walls, in the j- and ^-directions. The residual was reduced approximately 4.5 orders of 
magnitude over 5400 time steps. Navier execution times required for the three gas models were 
7.9-10 -5 , 1.1-10 -4 , and 1.1-10 -4 cpu-seconds per mesh-cell-point per time-step-iteration, 
respectively. 

Computed wall pressure distributions for turbulent comer flow are compared to 
experimental data [28] in Fig. 14. The present results again show good correlation with the data. 
Compared to the previous laminar comer flow predictions, these results indicate that the proper 
trending is produced by use of the two-comer wall turbulence model. 

CONCLUDING REMARKS 

The results presented herein show good correlation for all of the test cases considered. 
Since the equilibrium gas flux splitting schemes make use of the "equivalent" specific heat ratios, y 
and r, which are stored in the ^-vector as additional dependent variables, these results validate not 
only the implementation of the flux difference and flux vector splitting schemes, but also the 
restructuring of the present code to permit an arbitrary number of independent and dependent 
variables. 
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Both the Srinivasan and Tannehill equilibrium air model and the Liu and Vinokur 
generalized equilibrium gas model reproduce perfect gas results or, where appropriate, exhibit the 
proper real gas trends. With full vectorization, the equilibrium gas calculations were possible with 
only a small (-20%) increase in execution time. Successful coupling of the equilibrium 
air/equilibrium gas models with the one- or two-wall algebraic turbulence model was also 
demonstrated. 

Although not all of the code's capabilities were exercised, the results are indicative of the 
success of a substantial portion of the current effort. The resulting method should prove to be a 
valuable tool for use by the National Aero-Space Plane program, as well as a good starting point 
for future efforts aimed at incorporating nonequilibrium/chemistry effects, more sophisticated 
turbulence and transition models, or a variety of other physical phenomena. 
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(b) Temperature 

Figure 1: Supersonic Flat Plate Boundary Layer, Profiles at x = lm. 
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Figure 4: Hypersonic Flat Plate Boundary Layer, Surface Distributions. 
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(c) Pressure 

Figure 4: Hypersonic Flat Plate Boundary Layer, Surface Distributions (Conc’d). 
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(b) Pressure 

Figure 5: High Speed Inlet, Wall Distributions; Flux Difference Splitting. 
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X 

(c) Temperature 

Figure 5: High Speed Inlet, Wall Distributions; Flux Difference Splitting (Conc’d). 



X 


(a) Density 

Figure 6: High Speed Inlet, Wall Distributions; Flux Vector Splitting. 


25 






26 




HEAT TRANSFER, QDOT (W/m 2 x10~ 



Figure 7: Schematic of Computational Grid for Bent Nose Biconic. 



Figure 8: Bent Nose Biconic, Surface Heat Transfer Distribution. 
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Figure 9: Schematic of Computational Grid for Flared Cone; Laminar Flow. 



Figure 10: Flared Cone, Surface Distributions; Laminar Flow. 
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AXIAL LOCATION, X/L (L=24in) 

(b) Shear Stress 


12 r 

M oo = 1 6.93 Re=1. 976x1 0 5 /ft 


jCT 10 - CFL3D (perfect gas) O 

2 CFL3D (equilibrium air) 



AXIAL LOCATION. X/L (L=24in) 


(c) Pressure 

Figure 10: Flared Cone, Surface Distributions; Laminar Flow (Conc’d). 
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HEAT TRANSFER, QDOT (Btu/sec-ft 2 x1 CT 1 ) 


Figure 11: Schematic of Computational Grid for Flared Cone; Turbulent Flow. 



(a) Heat Transfer 


Figure 12: Flared Cone, Surface Distributions; Turbulent Flow. 




AXIAL LOCATION, X/L (L=24in) 

(b) Shear Stress 



AXIAL LOCATION. X/L (L=24in) 

(c) Pressure 

Figure 12: Flared Cone, Surface Distributions; Turbulent Flow (Conc’d). 
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distance from corner, (z— z woM )/x 

Figure 13: Laminar Corner Flow, Wall Pressure Distribution. 



distance from corner, (z-z W0 |,)/x 


Figure 14: Turbulent Corner Flow, Wall Pressure Distribution. 
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APPENDIX A 

Modified CFL3D Input Data File Description 
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*oom*«****m*om4m******* LINE TYPE ONE ****************************** 
title describing case 

*************************** LINE TYPE ONE. FIVE *************************** 
{DATA FOR LINE TYPE ONE. FIVE REPEATED FOR EACH FILE) 

path/name of binary grid file (unit 01) 

path/name of binary restart file (unit 02) 

path/name of binary PLOT 3D grid file (unit 03) 

path/name of binary PLOT3D flowfield file (unit 04) 

path/name of binary Liu & Vinokur equilibrium air coefficient file 

(unit 07) 

path/name of primary output file (unit 11) 

path/name of FIXI/FIXJ output file (unit 12) 

path/name of wing pressure output file (unit 14) 

path/name of secondary output file (unit 15) 

path/name of flowfield output file (unit 17) 

path/name of unsteady cp output file (unit 20) 

****************************** LINE TYPE TWO ****************************** 
xmach - freestream Mach number 

alpha - angle of attack 

beta - side-slip angle 

reue - freestream Reynolds number per unit length (millions) 

tinf - freestream temperature (degrees Rankine) 

isnd - wall temperature boundary condition flag 

- 0 adiabatic wall temperature 
• 1 specified wall temperature 

c2spe - wall temperature (temperature at wall divided by tmperature 

of freestream) 

if c2spe<-0, c28pe taken as freestream stagnation temperature) 

***************************** line TYPE THREE ***************************** 
sref - reference area 
cref - reference length 
bref - reference span 
xmc - moment center in x-direction 
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ymc - moment center in y-direction 
zmc - moment center in z-direction 


************************** LINE TYPE THREE. FIVE *************************** 

igas - perfect gas/equilibrium air flag 

- 1 perfect gas 

- 2 Tannehill equilibrium air 

- 3 Liu & Vinokur equilibrium air 

gamma - perfect gas, ratio of specific heats 
rgas - perfect gas, gas constant 
prgas - perfect gas, prandtl number 

scalex - meters per unit length 

- 1.0 if using meters 

- 0.3048 if using feet, default 

- 0.02540 if using inches 

scalet - degrees Kelvin per unit degree 

- 1.0 if using degrees Kelvin 

- 0.5556 if using degrees Rankine, default 

scalem - kilograms per unit mass 

- 1.0 if using kilograms 

- 14.59 if using slugs, default 

- 0.453472 if using pounds (mass) 


****************************** line type four ***************************** 

dt - time step 

< 0 local time stepping, CFL-abs(dt) 

> 0 constant time step (**dt) 

irest - 0 no restart 
- 1 restart 

iflagts - 0 constant dt 

>0 dt ramped over iflagts steps to dt*fmax 

fmax - maximum increase in dt 

iunst - 0 steady 

- 1 sinusoidal plunging 

- 2 sinusoidal pitching 

rfreq - reduced frequency 
alphau - Etching alpha 
cloc - pitching center 


****************************** LINE TYPE FIVE ***************************** 
ngrid - number of grids input 

nplot3d - number of flowfield data sets to be written in 
plot3d format 
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nprint - number of data seta to be sent to an output file 

nwrest - number of iterations between updates of the binarv 
restart file 1 


****************************** LINE TypjE SIX ********************„ ## * ##<k## 
{DATA FOR LINE TYPE SIX REPEATED NGRID TIMES) 


ncg - number of coarser grids to construct for multigrid/mesh 
sequencing {- 0 for embedded mesh) 


iem - embedded mesh flag 
- 0 for global grid 

■ 1 level of this embedded grid above global grid level 


iadvance - flag to skip any residual/update calculations 
>-0 proceed as usual 
< 0 skip residual/update calculations 


iforce - flag to skip the force routine 
>-0 proceed as usual 
< 0 skip force calculations 


imesh - mesh flag for grids topologically similar to: 

- 0 no singularties in mesh 

- 1 delta wing (AIAA 87-0207) 

- 2 prolate spheroid {AIAA 87-2627CP) 

- 3 prolate spheroid with sting (AIAA 87-2627CP) 

" 10 wing (o-h) 

* 11 wing (c-h) (AIAA 86-0274) 

- 12 wing (c-o) (AIAA 86-0274) 


ivisc (m) 


viscous/inviscid interaction flag m- 1 

0 inviscid 2 

1 laminar 3 

2 turbulent 


I-direction 

J-direction 

K-direction 


NOTE: The thin layer viscous terms can be included in either the 

f indirect ions f separately. The viscous terms can 

be included simultaneously in, at most, two directions, 
either j-k or i-k f for any particular grid. It is prefer- 
able to let k be the primary viscous direction and j be the 
secondary viscous direction. 


***************************** line TypE S£ven ***************************** 
(DATA FOR LINE TYPE SEVEN REPEATED NGRID TIMES) 


grid dimensions: 
idim - number pf 

* for imesh 

* for imeeh 

* for imesh 

* for imesh 

* for] jimesh 


points in i-direction 
* 1 axial direction (along chord) 

- 2 circumferentially along body 

“ 3 cir. along body/sting 

“ 10, 11 spanwise direction 

- 12 spanwise: wrapping around wing tip 


(h-mesh) 

(o-mesh) 

(c-mesh) 

(h-mesh) 

(o-mesh) 


jdim - number of points in j-direction 

* for imesh - 1,2,3 circumferentially along body/wing (c-mesh) 

* for imesh - 10 circumferentially along chord (o-mesh) 

* for imesh - 11,12 cir. along wing chord and wake (c-mesh) 

kdim - number of points in k-direction 

* for all imesh, radial direction 


itel - I location on body 
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- 

i at apex 




for 

imesh 


1 


- 

1 




for 

imesh 

- 

2, 3, 10,11, 12 

ite2 


i location on 

body 








- 

i at trailing 

edge 



for 

imesh 

- 

1 


- 

idim 




for 

imesh 

- 

2,3,12 


- 

i at wing tip 




for 

imesh 

- 

10,11 

jtel 

- 

j location on 

body 








- 

1 




for 

Imesh 

- 

1,2,3,10 


- 

j at trailing 

edge on 

lower 

surface 

for 

imesh 

- 

11,12 

jte2 

- 

j location on 

body 








- 

jdim 




for 

Imesh 

- 

1,2,3,10 


- 

j at trailing 

edge on 

upper 

surface 

for 

imesh 

- 

11,12 


**************************** LINE TYPE EIGHT ****************************** 
(DATA FOR LINE TYPE EIGHT REPEATED NGRID TIMES) 

inewg - restart flag for grid (not needed if ireet-0) 

- 0 read flowfield data from restart file 

- 1 initialize at freeatream or by linear interpolation 

from coarser grids 

igridc - grid to which this grid connects (input 0 for global 
mesh(iem-O) and the grid number in which the embedded 
mesh fits for embedded meshes (iem>0) ) 

js,ks,is - starting indices in connecting grid for placement of 
embedded mesh (input 0 for global meshes) 

je,ke,ie - ending indices in connecting grid for placement of 
embedded mesh (input 0 for global meshes) 

NOTE: The embedded meshes must be a regular refinement in all 

directions of the grid to which it connects. 


***************************** LINE TYPE NINE ****************************** 
(DATA FOR LINE TYPE NINE REPEATED NGRID TIMES) " 

idiag(m) - matrix inversion flag 

- 0 5x5 block tridiagonal inversion 

- 1 scalar tridiagonal inversions (recommended) 

iflim(m) - flux limiter flag m-1 : I-direction 

- 0 unlimited -2 : J-direction 

- 1 smooth limiter -3 : K-direction 

- 2 min-mod scheme (recommended) 


****************************** LINE TYPE TEN ****************************** 
] (DATA FOR LINE TYPE TEN REPEATED NGRID TIMES) 

ifds(m) - spatial differencing parameter for Euler fluxes 

- 0 flux-vector splitting 

- 1 flux-difference splitting (Roe's scheme) (recommended) 

rkapO (m) - spatial differencing parameter for Euler fluxes 

- -1 fully upwind 

- 0 Frommes's scheme 

- 1 central 

- 1/3 upwind-biased third order (recommended) 
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***************************** LINE T yp£ eleven **************************** 
(DATA FOR LINE TYPE ELEVEN REPEATED NGRID TIMES) 


boundary condition flags: 

mtypei(l) - boundary flag for i-0 boundary 
mtypei(2) - boundary flag for i-idim boundary 
mtypej(l) - boundary flag for j-0 boundary 
mtypej(2) - boundary flag for j-jdim boundary 
mtypek { 1 ) - boundary flag for k-0 boundary 
rotypek(2) - boundary flag for k — kdim boundary 


NOTE: Particular choices of mtypei/j/k determine the type of 

boundary conditions used at the edges of the computational 
grids and are best determined by inspection of subroutine 
BC. Additional boundary condition types can be 
incorporated into the algorithm by modifying subroutine BC 
according to the conventions outlined there* 


********************* ****** LINE TYPE ELEVEN ♦ ONE 
nbli - number of block boundary conditions 




*********** 


**************** LINE TypE eleven#two ************************** 
(DATA FOR LINE TYPE ELEVEN .TWO REPEATED NBLI TIMES) 


nblon - block boundary condition on or off ( >-0 or <0 ) 


************************** LINE TypE ELEVEN .three ************************* 

(DATA FOR LINE TYPE ELEVEN. THREE REPEATED NBLI TIMES) 


blckl 

ist 

jst 

kst 

ind 

jnd 

knd 

indl 

ind2 

blck.2 

ist 

jst 

kst 

ind 

jnd 

knd 

indl 

ind2 


- first block involved in block interface nbli 

- starting i-indice for blckl interface 

- starting j-indice for blckl interface 

- starting k-indice for blckl interface 

- ending i-indice for blckl interface 

- ending j-indice for blckl interface 

- ending k-indice for blckl interface 

- first indice which varies along blckl interface 
( 1-i ; 2-j ; 3-k ) 

- second indice which varies along blckl interface 
( 1-i ; 2-j ; 3-k ) 


second block involved in block interface nbli 
starting i-indice for blck2 interface 
starting j-indice for blck2 interface 
starting k-indice for blck2 interface 
ending i-indice for blck2 interface 
ending j-indice for blck2 interface 
ending k-indice for blck2 interface 
first indice which varies along blck2 interface 
< 1-i ; 2-j ; 3-k ) 


second indice which varies along blck2 interface 
( 1-i ; 2-j ; 3-k ) 




LINE TYPE TWELVE **************************** 


mseq - mesh sequencing flag for global grids (maximum 5) 

- 1 single solution on finest grid 

— 2 solution on second finest grid advanced ncyc(l) cycles 

followed by ncyc(2) cycles on finest grid. The solu- 
tion on the finest grid is obtained by interpolation 
from the coarser grid. If ncyc(2)-0, solution 
terminated on second finest grid after ncyc(l) steps 
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with restart file written for second finest grid at 
that point. 

> 2 sequencing from coarest to finest mesh as above 

mgflag - multigrid flag 

- 0 no multigrid 

- 1 multigrid on coarser global meshes 

- 2 multigrid on coarser global meshes and on 

embedded meshes 


iconsf - conservation flag 

- 0 nonconservative flux treatment for embedded grids 

- 1 conservative flux treatment for embedded grids 


mtt - 0 no additional iterations on the "up" portion 
of the multigrid cycle 

> 0 mtt additional iterations on the "up” portion 
of the multigrid cycle 


ngam - multigrid cycle flag 

- 1 V-cycle 

- 2 W-cycle 


*********** *****★*★**«★***** line TYPE THIRTEEN *************************** 
(REPEATED FOR EACH SEQUENCE 1 THROUGH MSEQ (COARSEST TO FINEST) ) 

ncycl - number of cycles 

mglevg - number of grids to use in multigrid cycling for 
the global meshes 

- 1 for single grid 

- 2 for two levels 

- m for m levels 

nemgl - number of embedded grid levels above the finest 

global grid (- 0 for global grids coarser than the 
finest global grid) 

- 0 no embedded grids 

- 1 one embedded grid 

-mm embedded grids 

nitfol - number of first order iterations 


A*************************** LINE TYPE FOURTEEN *************************** 
(REPEATED FOR EACH SEQUENCE 1 THROUGH MSEQ (COARSEST TO FINEST)) 

mitL - iterations on level L for each level L from coarsest 
to finest (mitL-1 recommended) 


ft*************************** LINE TYPE FIFTEEN **************************** 

(REPEATED NPL0T3D TIMES) 

block - designated block number for output 
istart - starting location in i-direction 
lend - erfding location in i-direction 
iinc - increment factor in i-direction 
jstart - smarting location in j-direction 
jend - ending location in j-direction 
jinc - increment factor in j-direction 
kstart - starting location in k-direction 
kend - ending location in k-direction 
kinc - increment factor in k-direction 


**************************** LINE TYPE SIXTEEN **************************** 

(REPEATED NPRINT TIMES) 
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block - 

istart 

lend 

line 

jstart 

jend 

jinc 

kstart 

kend 

kinc 




designated block number for output 

- starting location in i-direction 

- ending location In i-direction 

- increment factor in i-direction 

- starting location in j-direction 

- ending location in j-direction 

- increment factor in j-direction 

- starting location in k-direction 

- ending location in k-direction 

- increment factor in k-direction 
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Sample Modified CFL3D Input Data Files 
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Supersonic Laminar Flat Plate Boundary Layer 


2-d plate, cfl3dn - liu, lam 
binary grid file 

9 /scr6/rosen/plt32/plt32 .grd' 
binary restart file 

' /scr6/ rosen/plt 32/plt 32c .bin' 
plot3d binary grid file 

9 /scr 6/ rosen/plt 32 /pit 32c . pig' 
plot3d binary flowfield file 

9 /scr6/rosen/plt32/plt32c.plq' 

Liu 6 Vinokur binary equilibrium air coefficient file 
' / scr6/rosen/cf 13dn/ liu/ Hubs r . cof ' 
primary output file 

9 /scr 6 /rosen/plt 32 /pit 32c .out' 
fixi/fixj output file 

' /scr 6 /rosen/plt 32 /pit 32c . fix' 
wing pressure output file 

' /scr6/ rosen/plt 32 /pit 32c . wng' 
secondary output file 

' / scr 6 /rosen/plt 32 /pit 32c . sec' 
flowfield output file 

' /scr 6 /rosen/plt 32 /pit 32c.prt' 
unsteady cp output file 

' /scr 6 /rosen/plt 32 /pit 32c .ucp' 


XMACH 

ALPHA 

BETA 

REUE, MIL 

TINF, dK 

ISND 

C2SPE 


2.00 

0.000 

0.0 

1.650000 

221.60 

1 

1.0000 


SREF 

CREF 

BREF 

XMC 

YMC 

ZMC 



1.0000 

1.0000 

1.0000 

0 . 

0 . 

0 . 



I GAS 

GAMMA 

RGAS 

PRGAS 

SCALEX 

SCALET 

SCALEM 


3 

1.4 

286.9 

0.72 

1.0 

1.0 

1.0 


DT 

IREST 

IFLAGTS 

FMAX 

IUNST 

RFREQ 

ALPHAU 

CLOC 

- o.oox 

0 

500 

10.00 

0 

0.00000 

0.00000 

0.00000 

NGRID 

NPLOT3D 

NPRINT 

NWREST 





1 

0 

2 

100 





NCG 

IEM 

I ADVANCE 

I FORCE 

IMESH 

IVISCU) 

IVISC(J) 

IVISC(K) 

0 

0 

0 

0 

0 

0 

0 

1 

IDIM 

JDIM 

KDIM 

ITE1 

ITE2 

JTE1 

JTE2 


2 

51 

100 

1 

2 

1 

51 


INEWG 

IGRIDC 

IS 

JS 

KS 

IE 

JE 

KE 

1 

0 

0 

0 

0 

0 

0 

0 

IDIAG(I) 

IDIAG(J) 

IDIAG (K) 

IFLIM(I) 

IFLIM(J) 

IFLIM(K) 



0 

0 

0 

2 

2 

2 



IFDS (I) 

IFDS (J) 

IFDS <K) 

RKAP0 (I) 

RKAPO(J) 

RKAP0 (K) 



1 

1 

1 

0.33333 

0.33333 

0.33333 



MTYPEI (1) 

MTYPEI (2) 

MTYPEJ ( 1 ) 

MTYPEJ (2 ) 

MTYPEK (1) 

MTYPEK (2) 



11 

11 

27 

27 

67 

67 




NUMBER OF BLOCK INTERFACE BOUNDARY CONDITIONS 
0 


BLOCK INTERFACE BOUNDARY CONDITION ON OR OFF ( >-0 OR <0 ) 
BLCK1 — 


. 1ST 

JST KST IND 

JND KND 

IND1 IND2 

BLCK2 1ST 

MSEQ 

MG^LAG 

ICONSF 

MTT 

NGAM 

1 

0 

0 

0 

01 

NCYC 

MGLEVG 

NEMGL 

NITFO 


500 

01 

00 

000 


MIT1 

MIT2 

MIT3 

MIT4 

MIT5 

01 

01 

01 

01 

01 


PRINT OUT: 
BLOCK I START 


I END 
1 
1 


I INC JSTART 
1 50 

1 2 


JEND 

50 

50 


JINC KSTART 


KEND 

99 

1 


IND1 IND2 


KING 

1 

1 
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Hypersonic Laminar Flat Plate Boundary Layer 


2-d plate, cfl3dn : liu, lam 
binary grid file 

' /scr6/rosen/plt20/plt20 .grd' 
binary restart file 

' /scr6/rosen/plt20/plt20c .bin' 
plot3d binary grid file 

9 /scr6/rosen/plt 20 /pit 20c .pig' 
plot3d binary flowfield file 

9 /scr 6/rosen/plt 20 /pit 20 c .plq' 

Liu & Vinokur binary equilibrium air coefficient file 
9 /scr6/rosen/cf 13dn/liu/liubsr . cof ' 
primary output file 

'/scr6/rosen/plt20/plt20c. out' 
fixi/fixj output file 

'/scr6/rosen/plt20/plt20c.fix' 
wing pressure output file 

9 /scr6/rosen/plt20/plt20c . wng' 
secondary output file 

' /scr 6/rosen/plt 20/plt 20c . sec' 
flowfield output file 

r /scr6/rosen/plt20/plt20c.prt' 
unsteady cp output file 

' /scr 6/rosen/plt 20/plt 20c . ucp' 


XMACH 

ALPHA 

BETA 

REUE, MIL 

TINF, DR 

ISND 

C2SPE 


20.00 

0.000 

0.0 

0.200000 

100.00 

1 

10.0000 


SREF 

CREF 

BREF 

XMC 

YMC 

ZMC 



0.1000 

1.0000 

0.1000 

.05 

.5 

0 . 



IGAS 

GAMMA 

RGAS 

PRGAS 

SCALEX 

SCALET 

SCALEM 


3 

1.4 

286.9 

0.72 

1.0 

1.0 

1.0 


DT 

IREST 

IFLAGTS 

FMAX 

IUNST 

RFREQ 

ALPHAU 

CLOG 

-0.001 

0 

1500 

1000.00 

0 

0.00000 

0.00000 

0.00000 

NGRID 

NPLOT3D 

NPRINT 

NWREST 





1 

0 

2 

250 





NCG 

IEM 

I ADVANCE 

I FORCE 

IMESH 

IVISC(I) 

IVISC(J) 

IVISC(K) 

0 

0 

0 

0 

0 

0 

0 

1 

IDIM 

JDIM 

KDIM 

ITE1 

ITE2 

JTEl 

JTE2 


2 

65 

65 

1 

2 

1 

65 


INEWQ 

IGRIDC 

IS 

JS 

KS 

IE 

JE 

KE 

1 

0 

0 

0 

0 

0 

0 

0 

ID I AG ( I ) 

IDIAG (J) 

IDIAG (K) 

IFLIM(I) 

IFLIM(J) 

IFLIM(K) 



0 

0 

0 

2 

2 

2 



IFDS (I) 

IFDS (J) 

IFDS (K) 

RKAP0 { I ) 

RKAP0 (J) 

RKAP0 (K) 



1 

1 

1 

0.33333 

0.33333 

0.33333 



MTYPEI (1) 

MTYPEI (2) MTYPEJ <1 ) 

MTYPEJ (2) 

MTYPEK (1 ) 

MTYPEK (2) 



11 

11 

27 

27 

67 

67 



NUMBER OF 

n 

BLOCK INTERFACE BOUNDARY CONDITIONS 




BLOCK INTERFACE BOUNDARY CONDITION ON OR OFF ( >-0 

i OR <0 ) 



BLCK1 1ST 

JST KST IND JND KND 

IND1 IND2 

BLCK2 1ST 

JST KST IND JND KND 

IND1 IND2 

MSEQ 

MGFLAG 

ICONSF 

MTT 

NGAM 




1 

0 

0 

0 

01 




NCYC 

MGLEVG 

NEMGL 

NITFO 





1500 

01 

00 

000 





MIT 1 

MIT2 

MIT3 

MIT4 

MIT5 




01 

01 

01 

01 

01 




PRINT OUT: 







BLOCK I START IEND 

I INC JSTART JEND JINC KST ART KEND KINC 


1 

2 2 

1 

64 

64 1 

1 

65 2 


1 

2 2 

1 

2 

64 2 

1 

1 1 
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2-D INLET - AIAA 87-1117 (LIU, VAN LEER) 
binary grid file 

' /scr2/rosen/in/in . grd' 
binary restart file 

' /scr2/rosen/in/i3 .bin' 
binary plot3d grid file 
9 / scr2/rosen/in/ i3 .pig' 
binary plot3d flowfield file 
9 / scr2/rosen/in/i3 .plq' 

Liu 6 Vinokur binary equilibrium air coefficient file 
9 /scr 2/ rosen/cf 13dn/liu/liubsr , cof 9 
primary output file 

' /scr2/rosen/in/i3 .out' 
fixi/fixj output file 

9 /scr2/rosen/in/i3 . f ix # 
wing pressure output file 
' /scr2/rosen/in/i3 . wng' 
secondary output file 

9 /scr2/rosen/ in/ 13 . sec' 
flowfield output file 

' /scr2/rosen/in/i3 . prt ' 
unsteady cp output file 
' /scr2/rosen/in/i3 . ucp' 


XMACH 

ALPHA 

BETA 

REUE, MI L 

TINF, DK 

ISND 

C2SPE 


5.000 

0.000 

0.0 

4.940578 

3573.0 

1 

1.0 


SREF 

CREF 

BREF 

XMC 

YMC 

ZMC 



0.1 

1.0 

0.1 

0.5 

0.05 

0 . 



IGAS 

GAMMA 

RGAS 

PRGAS 

SCALEX 

SCALET 

SCALEM 


3 

1.4 

286.9 

0.72 

1.0 

1.0 

1.0 


DT 

IREST 

IFLAGTS 

FMAX 

IUNST 

RFREQ 

ALPHAU 

CLOC 

-0.010 

0 

300 

10.00 

0 

0.00000 

0.00000 

0.00000 

NGRID 

NPLOT3D 

NPRINT 

NWREST 





1 

0 

1 

100 





NCG 

IEM 

I AD VANCE 

I FORCE 

IMESH 

IVISC(I) 

IVISC(J) 

IVISC(K) 

0 

0 

0 

0 

1 

0 

0 

0 

IDIM 

JDIM 

KDIM 

ITE1 

ITE2 

JTE1 

JTE2 


2 

201 

51 

1 

2 

1 

201 


INEWG 

IGRIDC 

IS 

JS 

KS 

IE 

JE 

KE 

1 

0 

0 

0 

0 

0 

0 

0 

IDIAG(I) 

IDIAG(J) 

IDIAG(K) 

IFLIM(I) 

IFLIM(J) 

IFLIM(K) 



0 

0 

0 

2 

2 

2 



IFDS (I) 

IFDS (J) 

IFDS (K) 

RKAPO(I) 

RRAP0 (J) 

RKAPO(K) 



0 

0 

0 

0.33333 

0.33333 

0.33333 



MTYPEI (1) 

MTYPEI (2) 

MTYPEJ(l) 

MTYPEJ (2 ) 

MTYPEK (1 ) 

MTYPEK (2) 



27 

27 

27 

27 

27 

27 




NUMBER OF BLOCK INTERFACE BOUNDARY CONDITIONS 
0 


BLOCK INTERFACE BOUNDARY CONDITION ON OR OFF ( >-0 OR <0 ) 

BLCK1 1ST JST KST IND JND KND IND1 IND2 BLCK2 1ST JST KST IND JND KND IND1 IND2 


MSEQ 

MGf)I,AG 

ICONSF 

MTT 

NGAM 

1 

0 

0 

0 

01 

NCYC 

MGBEVG 

NEMGL 

NITFO 


300 

01 

00 

000 


MIT1 

MIT2 

MIT3 

MIT4 

MIT5 

01 

01 

01 

01 

01 


PRINT OUT: 

BLOCK ISTART I END I INC JSTART JEND JINC KSTART KEND KINC 

12211 201 1111 
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Bent Npg<? Bkonk 


BENT-BICONIC AT LOW-RE ALPHA-0 (NASA-TP-2334 ) 


binary grid file 

9 /scr/roaen/bnb/bnb . grd' 

binary restart file 

9 / gcr/ rosen/bnb/b3 . bin 9 
binary plot3d grid file 

' /9cr/rosen/bnb/b3 .pig' 

binary plot3d flowfield file 
' /scr /rosen/bnb/b3 . plq' 

Liu & Vinokur binary equilibrium air coefficient file 
9 liu/liubsr . cof ' 


primary output file 

9 / scr /rosen/bnb/b3 . out 9 
fixi/fixj output file 

9 /9cr/ro9en/bnb/b3 .fix' 

wing pressure output file 
9 / 9 cr /rosen/bnb/b3 ,wng' 
secondary output file 

9 /scr /rosen/bnb/b3 .sec' 
flowfield output file 


' /scr/rosen/bnb/b3 .prt' 
unsteady cp output file 


9 /scr/ rosen/bnb/b3 

.ucp' 






XMACH 

ALPHA 

BETA 

REUE, MIL 

TINF, DK 

ISND 

C2SPE 


9.860 

0.000 

0.0 

1.842000 

49.75 

1 

6.030151 


SREF 

CREF 

BREF 

XMC 

YMC 

ZMC 



0.001013 

0.121680 

0.121680 

0.067950 

0 . 

0 . 



IGAS 

GAMMA 

RGAS 

PRGAS 

SCALEX 

SCALET 

SCALEM 


3 

1.4 

286.9 

0.72 

1.0 

1.0 

1 . o 


DT 

IREST 

IFLAGTS 

FMAX 

IUNST 

RFREQ 

ALPHAU 

CLOC 

-0.001 

0 

400 

10.00 

0 

0.00000 

0.00000 

0.00000 

NGRID 

NPLOT3D 

NPRINT 

NWREST 





1 

0 

0 

100 





NCG 

IEM 

I ADVANCE 

I FORCE 

IMESH 

IVISC(I) 

IVISC(J) 

IVISC(K) 

1 

0 

0 

0 

3 

0 

0 

1 

IDIM 

JDIM 

KDIM 

ITE1 

ITE2 

JTE1 

JTE2 


85 

23 

45 

1 

85 

1 

23 


INEWG 

IGRIDC 

IS 

JS 

KS 

IE 

JE 

KE 

1 

0 

0 

0 

0 

0 

0 

0 

IDIAG (I) 

IDIAG (J) 

IDIAG (K) 

IFLIM(I) 

IFLIM(J) 

IFLIM(K) 



0 

0 

0 

2 

2 

2 



IFDS (I) 

IFDS (J) 

IFDS (K) 

RKAP0 (I) 

RKAP0 ( J) 

RKAP0 (K) 



0 

0 

0 

0.33333 

0.33333 

0.33333 



MTYPEI (1) 

MTYPEI (2) 

MTYPEJ(l) 

MTYPEJ (2 ) 

MTYPEK (1) 

MTYPEK (2) 



33 

67 

3 

3 

67 

77 



NUMBER OF 

BLOCK INTERFACE BOUNDARY CONDITIONS 




u 

BI.nrK INTERFACE BOUNDARY CONDITION ON OR OFF ( >-G 

\ OR <0 ) 



BLCK1 1ST 

JST KST IND JND KND 

IND1 IND2 

BLCK2 1ST 

JST KST IND JND KND 

IND1 IND2 

MSEQ 

MGFLAG 

ICONSF 

MTT 

NGAM 




2 

0 

0 

0 

01 




NCYC 

MGLEVG 

NEMGL 

NITFO 





400 

01 

00 

000 





000 

01 

00 

000 





MIT 1 

MIT2 

MIT3 

MIT4 

MIT5 




01 

01 

01 

01 

01 




01 

01 

01 

01 

01 





PRINT 

BLOCK 

OUT: 

I START 

IEND 

I INC 

JSTART 

1 

1 

85 

1 

1 


JEND JINC KSTART KEND KINC 
23 111 1 
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A50-FLARED-CONE AT LOW-RE RUN-17 (AFFDL-TR-65-199) 
binary grid file 

' /scr2/rosen/a50/a50 . grd' 
binary restart file 

f /scr2/rosen/a50/a3.bin' 
binary plot3d grid file 

9 /scr2/rosen/a50/a3 . pig' 
binary plot3d flowfield file 
9 /scr2/rosen/a50/a3 .plq' 

Liu & Vinokur binary equilibrium air coefficient file 
' /scr2/rosen/cf 13dn/liu/liubsr . cof ' 
primary output file 

9 /scr2/rosen/a50/a3 . out' 
fixi/fixj output file 

' / sc r 2 /rosen/a 50 /a 3. fix' 
wing pressure output file 
' /scr2/rosen/a50/a3.wng' 
secondary output file 

9 /scr2/rosen/a50/a3 . sec' 
flowfield output file 

' /scr2/rosen/a50/a3 .prt' 
unsteady cp output file 

9 / sc r2 /rosen/a 50 /a3 . ucp' 


XMACH 

ALPHA 

BETA 

REUE, MIL 

TINF, DR 

ISND 

C2SPE 


16.930 

0.000 

0.0 

0.197600 

83.73 

1 

6.329870 


SREF 

CREF 

BREF 

XMC 

YMC 

ZMC 



1.0 

1.0 

1.0 

1.0 

0 . 

0 . 



IGAS 

GAMMA 

RGAS 

PRGAS 

SCALEX 

SCALET 

SCALEM 


3 

1.4 

1715.6 

0.72 

0.0 

0.0 

0.0 


DT 

IREST 

IFLAGTS 

FMAX 

IUNST 

RFREQ 

ALPHAU 

CLOC 

-0.001 

0 

300 

10.00 

0 

0.00000 

0,00000 

0.00000 

NGRID 

NPLOT3D 

NPRINT 

NWREST 





1 

0 

0 

100 





NCG 

IEM 

IADVANCE 

I FORCE 

IMESH 

IVISC(I) 

IVISC(J) 

IVISC(K) 

1 

0 

0 

0 

1 

0 

0 

1 

IDIM 

JDIM 

KDIM 

ITE1 

ITE2 

JTE1 

JTE2 


97 

19 

45 

5 

97 

1 

19 


INEWG 

IGRIDC 

IS 

JS 

KS 

IE 

JE 

KE 

1 

0 

0 

0 

0 

0 

0 

0 

IDIAG(I) 

IDIAG(J) 

IDIAG(K) 

IFLIM(I) 

IFLIM(J) 

IFLIM(K) 



0 

0 

0 

2 

2 

2 



IFDS { I ) 

IFDS ( J) 

IFDS (K) 

RKAPO(I) 

RKAPO(J) 

RKAP0 (K) 



1 

1 

1 

0.33333 

0.33333 

0.33333 



MTYPEI (1) 

MTYPEI (2) 

MTYPEJ { 1 ) 

MTYPEJ (2 ) 

MTYPEK(l) 

MTYPEK (2) 



67 

67 

1 

1 

67 

77 




NUMBER OF BLOCK INTERFACE BOUNDARY CONDITIONS 
0 

BLOCK INTERFACE BOUNDARY CONDITION ON OR OFF < >-0 OR <0 ) 


BLCK1 1ST JST KST IND JND KND INDl IND2 BLCK2 1ST JST KST IND JND KND IND1 IND2 


MSEQ 

MPFLAG 

ICONSF 

MTT 

NGAM 

2 

0 

0 

0 

01 

NCYC 

MGLEVG 

NEMGL 

NITFO 


300 

01 

00 

000 


000 

01 

00 

000 


MITl 

MIT2 

MIT3 

MIT4 

MIT5 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 


PRINT OUT: 

BLOCK I START IEND I INC JSTART JEND JINC KSTART KEND KINC 

11 97 11 19 1111 
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Flared Cone (Turbulent) 

A50- FLA RED-CONE AT HIGH-RE RUN-32 (AFFDL-TR-65-199) 
binary grid file 

' /acr/ro9en/a32/a32 . grd' 
binary restart file 

'/scr/rosen/a32/a3.bin' 
binary plot3d grid file 
'/9cr/rosen/a32/a3.plg' 
binary plot3d flowfield file 
f /scr/rosen/a32/a3 . plq' 

Liu & Vinokur binary equilibrium air coefficient file 
' liu/liubsr . cof ' 
primary output file 

f /scr/rosen/a32/a3 .out' 
fixi/fixj output file 

'/scr/rosen/a32/a3.fix' 
wing pressure output file 
9 /scr/ rosen/a32/a3 . wng' 
secondary output file 

'/scr/rosen/a32/a3.sec' 
flowfield output file 

'/scr/rosen/a32/a3 . prt' 
unsteady cp output file 
'/scr/ro9en/a32/a3.ucp' 


XMACH 

ALPHA 

BETA 

REUE, MIL 

TINF,DR 

ISND 

C2SPE 


7.850 

0.000 

0.0 

4.697000 

130.2 

1 

4.070661 


SREF 

CREF 

BREF 

XMC 

YMC 

ZMC 



1.0 

1.0 

1.0 

1.0 

0 . 

0 . 



IGAS 

GAMMA 

RGAS 

PRGAS 

SCALEX 

SCALET 

SCALEM 


3 

1.4 

1715.6 

0.72 

0.0 

0.0 

0.0 


DT 

IREST 

IFLAGTS 

FMAX 

IUNST 

RFREQ 

ALPHAU 

CLOC 

-0.010 

0 

300 

10.00 

0 

0.00000 

0.00000 

0.00000 

NGRID 

NPLOT3D 

NPRINT 

NWREST 





2 

0 

0 

100 





NCG 

IEM 

I ADVANCE 

I FORCE 

IMESH 

IVISC(I) 

IVISC(J) 

IVISC(K) 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

1 

IDIM 

JDIM 

KDIM 

ITE1 

ITE2 

JTE1 

JTE2 


5 

19 

45 

5 

5 

1 

19 


93 

19 

45 

1 

93 

1 

19 


INEWG 

IGRIDC 

IS 

JS 

KS 

IE 

JE 

KE 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

IDIAG(I) 

ID I AG ( J) 

I D I AG < K ) 

IFLIM(I) 

IFLIM(J) 

IFLIM(K) 



0 

0 

0 

2 

2 

2 



0 

0 

0 

2 

2 

2 



IFDS (I) 

IFDS (J) 

IFDS (K) 

RKAP0 (I) 

RKAP0 (J) 

RKAP0 (K) 



1 

1 

1 

0.33333 

0.33333 

0.33333 



1 

1 

1 

0.33333 

0.33333 

0.33333 



MTYPEI (1) 

MTYPEI (2) 

MTYPEJ(l) 

MTYPEJ (2 ) 

MTYPEK (1) 

MTYPEK (2) 



67 

67 

1 

1 

67 

77 



67 

67 

1 

1 

67 

77 




NUMBER OF BLOCK INTERFACE BOUNDARY CONDITIONS 


1 

BLOCK INTERFACE BOUNDARY CONDITION ON OR OFF ( >-0 OR <0 ) 

1 

BLCK1 1ST JST *KST IND JND KND IND1 IND2 BLCK2 1ST JST KST IND JND KND IND1 IND2 


5 

115 

19 45 

2 3 

2 1 

MSEQ 

MGFLAG 

ICONSF 

MTT 

NGAM 

2 

0 

0 

0 

01 

NCYC 

MGLEVG 

NEMGL 

NITFO 


300 

01 

00 

000 


000 

01 

00 

000 


MIT1 

MIT2 

MIT3 

MIT4 

MIT5 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 
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Laminar Corner Flow 


symmetric wedge corner : liu,lam 
binary grid file 

' /scr2/rosen/corner/lam.grd' 
binary restart file 

' /scr2/rosen /corner/ 13 .bin' 
plot3d binary grid file 

' /scr 2 /rosen/corner/ 13 . pig' 
plot3d binary flowfield file 
9 /scr 2 /rosen/corner/ 13 .plq r 
Liu £ Vinokur binary equilibrium air coefficient file 
' liu/liubsr . cof ' 
primary output file 

' /scr 2/ rosen/corner/ 13 .out' 
fixi/fixj output file 

9 / scr 2 /rosen/corner/ 13 .fix' 
wing pressure output file 

' / scr2/rosen/corner/13 . wng' 
secondary output file 

' / scr 2 /r os en /corner/ 13 .sec' 
flowfield output file 

9 /scr 2 /rosen/corner/ 13 .prt ' 
unsteady cp output file 

' / scr2/rosen/corner/13 . ucp r 


XMACH 

ALPHA 

BETA 

REUE, MIL 

TINF, DK 

ISND 

C2SPE 


3.00 

0 . 000 

0.0 

3.07 

105.0 

1 

2.8 


SREF 

CREF 

BREF 

XMC 

YMC 

ZMC 



1.0000 

1.0000 

1.0000 

0.0 

0.0 

0.0 



IGAS 

GAMMA 

RGAS 

PRGAS 

SCALEX 

SCALET 

SCALEM 


3 

1.4 

286.9 

0.72 

1.0 

1.0 

1.0 


DT 

IREST 

IFLAGTS 

FMAX 

IUNST 

RFREQ 

ALPHAU 

CLOC 

-0.0X0 

0 

300 

10.00 

0 

0.00000 

0.00000 

0.00000 

NGRID 

NP LOT 3D 

NPRINT 

NWREST 





1 

0 

0 

300 





NCG 

IEM 

I AD VANCE 

I FORCE 

IMESH 

IVISC(I) 

IVISC(J) 

IVISC(K) 

0 

0 

0 

0 

0 

0 

1 

1 

IDIM 

JDIM 

KDIM 

ITE1 

ITE2 

JTE1 

JTE2 


2 

121 

121 

1 

2 

1 

121 


INEWG 

IGRIDC 

IS 

JS 

KS 

IE 

JE 

KE 

1 

0 

0 

0 

0 

0 

0 

0 

IDIAG(I) 

IDIAG(J) 

IDIAG(K) 

IFLIM(I) 

IFLIM(J) 

IFLIM(K) 



0 

0 

0 

2 

2 

2 



IFDS (I) 

IFDS (J) 

IFDS (K) 

RKAP0 (I) 

RKAPO(J) 

RKAP0 (K) 



1 

1 

1 

0.33333 

0.33333 

0.33333 



MTYPEI (1) 

MTYPEI (2) 

MTYPEJ ( 1 ) 

MTYPEJ (2) 

MTYPEK(l) 

MTYPEK (2) 



1002 

1002 

1004 

1002 

1004 

1002 




NUMBER OF BLOCK INTERFACE BOUNDARY CONDITIONS 
0 


BLOCK INTERFACE BOUNDARY CONDITION ON OR OFF ( >-0 OR <0 ) 

BLCK1 1ST JST KST IND JND KND IND1 IND2 BLCK2 1ST JST KST IND JND KND IND1 IND2 


MSEQ 

MGFLAG 

I CONS F 

MTT 

NGAM 

1 

0 

0 

0 

01 

NCYC 

MGLEVG 

NEMGL 

NITFO 


300 

01 

00 

000 


MIT1 

MIT2 

MIT3 

MIT4 

MIT5 

01 

01 

01 

01 

01 
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Turbulent Corner Flow 


symmetric wedge corner : liu r turbulent 
binary grid file 

' / scr2/rosen/corner/turb3 . grd' 
binary restart file 

'/ scr2/rosen/corner/t3.bin' 
plot3d binary grid file 

' /scr 2/rosen/cor ner/t 3 .pig' 
plot3d binary flowfield file 
' / scr2/rosen/ corner/t3 . plq' 

Liu & Vinokur binary equilibrium air coefficient file 
' liu/liubsr.cof' 
primary output file 

' / scr2/rosen/corner/t3 .out f 
fixi/fixj output file 

' /scr 2 /rosen/cor ner/t 3 .fix' 
wing pressure output file 

' /scr 2 /rosen/cor ner/t 3 .wng' 
secondary output file 

'/scr2/rosen/corner/t3.sec' 
flowfield output file 

'/scr2/rosen/corner/t3.prt' 
unsteady cp output file 

' / scr2/rosen/corner/t 3 . ucp' 


XMACH 

ALPHA 

BETA 

REUE, MIL 

TINF, DK 

ISND 

3.00 

0.000 

0.0 

3.2189 

105.0 

1 

SREF 

CREF 

BREF 

XMC 

YMC 

ZMC 

1.0000 

1.0000 

1.0000 

0.0 

0.0 

0.0 

IGAS 

GAMMA 

RGAS 

PRGAS 

SCALEX 

SCALET 

3 

1 . 4 

286.9 

0,72 

1.0 

1.0 

DT 

IREST 

IFLAGTS 

FMAX 

IUNST 

RFREQ 

-0.010 

NGRID 

1 

0 

NPLOT3D 

0 

300 

NPRINT 

0 

10.00 

NWREST 

300 

0 

0.00000 

NCG 

I EM 

I ADVANCE 

I FORCE 

IMESH 

IVISC(I) 

0 

0 

0 

0 

0 

0 

IDIM 

JDIM 

KDIM 

ITEl 

ITE2 

JTE1 

2 

121 

121 

1 

2 

1 

INEWG 

IGR1DC 

IS 

JS 

KS 

IE 

1 

0 

0 

0 

0 

0 

IDIAG (I) 

IDIAG ( J) 

IDIAG (K) 

IFLIM(I) 

IFLIM(J) 

IFLIM(K) 

0 

0 

0 

2 

2 

2 

IFDS (I) 

IFDS (J) 

IFDS (K) 

RKAP0 (I) 

RKAP0 (J) 

RKAP0 (K) 

1 

1 

1 

0.33333 

0.33333 

0.33333 

MTYPEI (1) 

MTYPEI (2) 

MTYPEJ(l) 

MTYPEJ (2 ) 

MTYPEK(l) 

MTYPEK (2 ) 

1002 

1002 

1004 

1002 

1004 

1002 

NUMBER OF 
n 

BLOCK INTERFACE BOUNDARY CONDITIONS 

BLOCK INTERFACE BOUNDARY CONDITION ON OR 

L OFF ( >-0 

OR <0 ) 

JST KST IND 

BLCK1 1ST 

JST rfST IND JND KND 

IND1 IND2 

BLCK2 1ST 

MSEQ 

MGFLAG 

I CONS F 

MTT 

NGAM 

1 

NCYC 

300 

0 

MGLEVG 

01 

0 

NEMGL 

00 

0 

NITFO 

000 

01 


MIT1 

MIT2 

MIT3 

MIT4 

MIT5 


01 

01 

01 

01 

01 



C2SPE 

2.8 


SCALEM 

1.0 

ALPHAU 

0.00000 


IVISC(J) 

1 

JTE2 

121 

JE 

0 


JND KND 


CLOC 

0.00000 


IVISC(K) 

1 


KE 

0 


IND1 IND2 
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